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Abstract 

We consider the contributions of both color-singlet and 
octet cc bound states for studying the suppression of 
J/ip production in quark gluon plasma (QGP) and have 
found that the existence of the color-octet cc may cause 
the pr— distribution oi J/ip a decrease with px after a 
critical value {pT)max, which is quite different from the 
behavior in Hadronic Matter (HM). 
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I. Introduction 

It was expected that a suppression of J /if) production in relativistic heavy ion 
collisions can serve as a clear signature for the formation of a new matter phase 
Quark Gluon Plasma (QGP) This suppression effect was observed by NA38 col- 
laboration later 0. However, successive research pointed out that such suppression 
could also exist in Hadronic Matter (HM), even though by a completely different 
mechanism 0] . Therefore it is believed that an observation of mere J / if) production 
rate suppression is not sufficient to confirm QGP formation, and a more thorough 
study on the characteristics of the J /if) suppression is necessary and may lead to 
final confirmation of QGP state. The pt distribution of J/ip may be a good candi- 
date, since the different mechanisms which result in J /if) suppression in both HM 
and QGP may provide different pt distributions of J/-?/' 0] . Our numerical 
results confirm this conjecture (see below). 

In this work, we discuss the px distribution of J/ip in QGP and analyze differ- 
ences from that in HM. 

In principle, the J/ip quarkonium is described in a Fock state decomposition 

\J/ij > = 0{l)\ccCs['^) > +0{v)\ccCpf)g > +0{v')\ccCsl>% > + 

0{v')\ccCs['^%g > +Oiv')\cci'D^j^%g > +..., (1) 

where 2^+1 ^(i>8) charaterizes the quantum state of the cc with color-singlet or octet 
respectively. This expression is valid for the non-relativistic QCD (NRQCD) frame- 
work and the coefficients of each components depend on the three- velocity |?T| of the 
heavy quark, and under the limit of l^l 0, i.e. both c and c remain at rest, eq.(|l]) 
recovers the expression for color-singlet picture of J/ip where 0(1) = 1. 

The color-octet component {cc)s has been proved to play an important role 
for interpreting the CDF experimental data p|,(see next section for more detail), 
therefore one can expect that the color-octet {cc)s would also manifest itself in heavy 
ion collisions. Since the color-octet can interact with gluons much more strongly 
than the color-singlet (cc)i, it would dissolve much faster into D and D than (cc)i, 
while traversing through a gluon-rich environment such as QGP, Moreover, the linear 
confinement potential ar at finite temperature turns out to be in the form 
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a{T)r = ao{^ )^0(T,,, - T)r, (2) 

dec 

where 5 is a parameter. So after T > Tdec, the hnear confinement disappears and 
the Casimir operator < Sl^^^lS >= +3 > 0, the Coulomb-type potential would 
become repulsive for an color-octet {cc)s- It means that in QGP c and c cannot be 
bound together in an octet, and the original color-octet bound state {cc)s which is 
produced at the early stage of the collision would dissolve. 

In other words, in principle, the color-octet {cc)s does exist in HM, but could 
not be built up in QGP. However, in our phenomenological study, the main fraction 
of {cc)s dissolves in QGP, while a small amount remains (see next section). 

In this scenario, the {cc)s does not make a substantial contribution in QGP, as it 
does in pp— collisions where the color-octet {cc)s constitutes the dominant part for 
large pt distribution of J/iJ. Thus one would expect that the large px contribution 
of J/i/j would be more suppressed than the smaller pt one in QGP, but this does 
not occur for HM and we will discuss it in more details later. 

In next section, we describe the color-octet scenario borrowed from the pp colli- 
sions and discuss a simple but popular model for J/ip suppression in QGP. In the 
third section we present the numerical results, while the last section is devoted to 
our conclusion and discussion. 

II. Formulation 

(i) The color-octet J/ip. 

It has already been demonstrated that the Fock decomposition form of J/ ip can 
be expressed as eq.(^ and the coefficients of each term are functions of the velocities. 
Calculating the coefficients has been proved to be beyond the present ability, since 
this is non-perturbative QCD process. In principle, by the information gained in pp 
collisions, the produced cc can be in either color-singlet or color-octet. According 
to the newly developed factorization formalism the inclusive J/ip production 
amplitude from the collision of partons a, b can be factorized into short- and long- 
distance pieces: 

A{ab ^ J/i) + X)= J2 Aab cc{^'^^Lf)X) x A{cc{^'^^Lf) ^ J/i;). (3) 

(c)L,J 
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The short-distance part for the hard production of cc pair can be calculated 
within perturbative QCD. On the other hand, the long-distance component of the 
hadronization cc(^^~^^L^j^) — > J/ip is a non-perturbative process, which can be cal- 
culated by lattice simulation or expressed in terms of some phenomenological pa- 
rameters. 

For the pt distribution, the differential cross section reads [0,|r^ 

da{pp^ J/lp + X) ^ da{pp {J/i))octet + X) ^ da{pp {J / 1p) singlet + X) 

dpT dpT dpT 

where the subscripts denote different color components of eq.(|I]). At the large pt 
limit, the main contribution comes from gluon fragmenting into a color-octet cc pair. 

It is important to note that in hadron collisions, the non-perturbative process of 
evolution from the produced cc to J/ip is expressed in terms of some phenomeno- 
logical parameters which are independent of px-, whereas in hadron- nucleus and 
nucleus- nucleus collisions, the hadronization of cc depends on px- In fact, we take 
the parameters achieved in hadron collision to be given and consider the dependence 
can be obtained automatically for ion collision situation(see below). 

With the understanding, we can assume that the color-singlet and octet cc are 
produced at early stage of heavy ion collisions and later evolve according to the 
environment. If the temperature and density have not reached the phase transition, 
the system remains in HM, but if the phase transition indeed occurs, a QGP phase 
is formed, thus the (cc)i and (cc)8 would have different evolution rules. It is noted 
that here we assume, the cc are produced at the early stage of the ion collisions 
through hard N-N collisions when the phase transition has not occurred yet, and 
then the system would either remains in HM or a new phase QGP is formed at later 
time. Thus we only need to deal with the evolution of (cc). Namely we ignore the 
production of cc by the deconfined quarks or gluons in QGP. 

(ii) A typical mechanism for suppression of J/ip production in QGP. 

The Pt distribution of J / ip can be described by a fundamental picture |^ |^ . Due 
to the Debye screening, the potential between c and c would undergo a modification 
from Vcoui = into ^e~'^^^° where aeff oc < 1, SjA'^A"!!, 8 >, and when vd is 
sufficiently small, (cc)i can no longer form J/ip. In addition, as discussed above, due 
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to the repulsive interaction between cc in octet state in QGP, the {cc)s produced at 
the early stage of the collision can in principle, never substantialize into J/tp in the 
QGP background if there is enough time or the size of the QGP region is sufficiently 
large. But we will see, it is not always the case. As for the case in HM case, (cc)i and 
{cc)s also have possibility to dissolve into DD pairs by scattering with the hadronic 
matter to reduce the J/ip production rate. It is noted that the mechanisms in HM 
and QGP for the dissolution are different, for the color-singlet and octet they lead 
to even more apart results. 

To concretely calculate the changes, let us take a reasonable and simplified phys- 
ical picture p . Provided that in A-A collisions, the initially produced bound 
state (cc)i has a very small radius and needs approximately ~ 0.7 fm to reach the 
intrinsic size of J/ tp. Thus the (cc)i bound state of small volume does not suffer from 
the Debye screening in QGP. In the laboratory frame, due to time dilation, one has 
t^p = 'Ti/.y^l + {px/My, which turns out to be the formation time assuming pz = 0. 
For the case of Pz 7^ 0, this result applies to the transverse distribution. Therefore 
if the bound state (cc)i emerges from the QGP region, within a time interval tf less 
than t^, it can ignore the QGP effect and eventually forms J/ip meson which can 
be observed experimentally. There are two models, the static and the plasma 
expansion [Q. Because the QGP phase is produced at RHIC and evolves very fast, 
we take the second model to study the pr distribution of J/ip. For {cc)s, as the 
potental is repulsive in QGP, {cc)s suffers the effect of QGP almost immidiently, 
so the corresponding time being the time for {cc)s starting to suffer from QGP 
effect should be much shorter than due to the strong interaction with the gluons 
in environment. 

Introducing Fi as the ratio of the number of produced (cc)i, which eventually 
forms J/tp to that of the total (cc)i produced in collisions. Similarly, we could also 
define a ratio Fg for {cc)s- Obviously 

Fi=0 tf tf>t^, 

(5) 

Fi = 1 if tf < t^, 

I. 

where t/ is the time as (cc)i lingers in QGP. By the Bjorken hydrodynamic model. 
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CdrJ,^ d<Prp{r)9m- tf{r, 
TT /q^° drrp{r) 



FiiPT) = ^ :,n:L\L (6) 



where p{r) = (1 — r"^ / Rq)''6{Rq — r) with b=l is the initial distribution density of 
(cc)i which is located at r from the center. It is easy to find that there is a maximum 
transervse momentum {pT)max and one has Fi = 1 as long as px > {PT)max- 



{pT)ma. = MJi — y-l, (7) 

V 

where tm is the lifetime of the QGP and M is the mass of the J/?/)— meson. It is 
noted that Fi = 1 as long as pr > {pT)max- Actually Eq.(7) only holds in cases 
where the plasma lives for a time short compared to the average transit time ti of 
a cc in the system. This is because of the finite size effect: a (cc)i with a large pt 
escape the interaction region before it is formed 0]. As has been discussed in Ref. 

|jl2|, in general, for A — A collisions the transverse expansion is comparatively 
small, so that the condition tm < ti is satisfied. 

When d = V ■ <^ R{t^) where v is the velocity of cc, one can obtain an 
approximate expression 



with approximately a = 1/2, b = 1 for equal nuclear collision. In this expression, 
the surface effects are neglected compared to eq.®. For not too large pt and large 
nuclei, the results by the two expressions are very close |^ Therefore, we will 
use this approximate one to make a qualitative discussion on the problem below. 



(iii) The effects of (cc)8 on the pt distribution. 

As understood, at the initial stage of the A-A collisions, just in analog to the 
hadron-hadron collisions, there are both (cc)i and (cc)8 produced, by the information 
gained from hadron collisions, their relative fractions are pT~dependent and can be 
expressed as fiipx) and fsipr) respectively (a simplified notation from eq.(P). They 
would also make different contributions to the pt distribution of J/ip- Thus the ratio 
of the number of (cc) which finally become J/ to the total number of originally 
produced (cc) could be expressed as 
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F = FMMpt) + Fs{PT)fs{pT) 



(9) 



where Fi and Fg are the corresponding ratios for (cc)i and {cc)s respectively as 
defined above. 

We have argued that since {cc)s is colored, it interacts strongly with gluons in the 
environment of QGP and can dissolve much faster than (cc)i which is color-blind. 
One can expect that the {cc)s does not contribute to px distribution much, namely 
Fg ^ Fi in QGP. Moreover, if the potential picture given in eq.(|l]) is valid, {cc)s 
cannot survive in QGP after all, thus Fg = 0, if the time the {cc)s needs to traverse 
cross the QGP region is very long, there would be no {cc)s contribution to J/ip. 
However, this has only a statistical meaning. For a finite size of the region, there 
should be fiuctuation that a small fraction of {cc)s emerging out. Anyhow, all models 
we have so far are phenomenological and are not derived from the first principle of 
QCD, therefore one can expect some small declination from that determined by the 
model. So even though unable to derive Fg, we can phenomenologically choose some 
reasonable ratios as Fg = Fi/2, Fi/10 and Fg = for our numerical computations. 

III. The numerical results 

With our simple model and phenomenological choice of Fg, we have calculated 
the pt distribution of J/ip. The numerical results are depicted in Fig.l. 

Because the {pT)max (see eq.(^) depends on the formation time of J/ip and the 
lifetime of the QGP, its estimation cannot be accurate in our present knowledge, 
so that we take it as a free parameter and the only one in the whole calculation. 
Thus one can notice that the position of {pT)max is fixed just because we take it as 
a parameter and everything is scaled with this value. We have found that changes 
of {pT)max do uot affcct our qualitative conclusion that we predict an increase of 
F with increase of pt at pt < {pT)max, but a linear decrease with further increase 
of Pt after px > {pT)max- The choices of the ration Fg/Fi also do not affect our 
qualitative conclusion. 

More discussions will be given in next section. 

IV. Conclusion and discussion 
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We employ a simple model to estimate the px distribution of J/ ip which emerges 
from a QGP region. The data of hadron collisions indicate that the {cc)s component 
is not negligible and especially plays a crucial role for the large pr distribution of 
produced J/ip. It motivates us to consider effects of the {cc)s component on the pt 
distribution of J/tp in relativistic heavy ion collisions, especially when the excited 
region turns into the QGP phase. Indeed, {cc)s produced at the early stage of 
collision would have a completely different behavior of evolution in QGP and HM. 
It is easy to understand because QGP phase provides a gluon-rich environment which 
strongly interacts with the color-octet {cc)s while there are no free gluons in the HM 
background, thus the interaction is much weaker. Therefore the pt distribution may 
be taken as a signature of QGP formation. 

Since the mechanism for J/ip suppression in QGP is obviously simpler than that 
in HM, so we first study the pt distribution of J/ ip emerging from the QGP region 
by calculating the F-value while only keep a qualitative analysis of corresponding 
observation for HM. 

Our numerical results and qualitative arguments all indicate that as px > {pT)max 
-F— value decreases with px and it is in contrary to the behavior of F at HM where 
F always increases with px- That is significant. It is believed that in the previous 
O — U and S — U experiments, the energy scale is not high enough to produce a 
QGP state, and then the observed F— value indeed keeps increasing with an increase 



of Pt- The hadron picture gives a satisfactory description to the result 1T2| [Q, so 
it confirms that a turning over of the dependence of F on pt may signify a QGP 
formation. 

An unusual suppression of J/ ip has been observed at present Ph — Pb collision 
experiment and there have been some discussions about the possible mechanisms 
fl^ . Based on our model, further measurement of the pr-dependence of the J/ip 
suppression may provide additional information about the formation of QGP. 

As discussed above, the suppression of J/ ip production also occurs at HM, so the 
phenomenon is hard to be taken as a clear signature of QGP, however, a detailed 
analysis on its characteristics may do the job. The produced (cc)i and {cc)s evolve 
according to different ways in HM and QGP, so that the pr distribution of J/ip 
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determined by their contributions would be different in two matter pliases. One can 
optimistically liope tliat an analysis of tlie px distribution of J/ip, especially the 
dependence of F— value on pt may inform us of existence of QGP. 

However, to make a convincing conclusion, one should do a quantitative estima- 



tion of F— value at HM and QGP and this is our goal of next work ||T5[. We believe 
that we cannot confirm the QGP formation even though we observe an F-decrease 
with Pt unless we have a thorough study on its behavior at HM. All these not only 
require more theoretical work which we will publish elsewhere, but also demand 
more accurate data which will be given by RHIC experiments. 
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Figure Captions 

Fig.l, the dependence of the F— values on px, where different curves correspond 
to various ratios of Fg/Fi which are phenomenologically introduced. The (pr)maa: is 
the only parameter in the scenario which is discussed in context. 
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